The discretised harmonic oscillator: Mathieu functions and a new class of generalised 

Hermite polynomials 



M. AuhoIeEI 

Department of Physics, University of Jyvdskyld, P.O. Box 35 (YFL), FIN-4OOI4 University of Jyvdskyld, FINLAND 

(Dated: February 7, 2008) 

We present a general, asymptotical solution for the discretised harmonic oscillator. The corre- 
sponding Schrodinger equation is canonically conjugate to the Mathieu differential equation, the 
Schrodinger equation of the quantum pendulum. Thus, in addition to giving an explicit solution 
for the Hamiltonian of an isolated Josephon junction or a superconducting single-electron transistor 
(SSET), we obtain an asymptotical representation of Mathieu functions. We solve the discretised 
harmonic oscillator by transforming the infinite-dimensional matrix-eigenvalue problem into an in- 
. . . finite set of algebraic equations which are later shown to be satisfied by the obtained solution. The 

, proposed ansatz defines a new class of generalised Hermite polynomials which are explicit functions 

■ of the coupling parameter and tend to ordinary Hermite polynomials in the limit of vanishing cou- 

' pling constant. The polynomials become orthogonal as parts of the eigenvectors of a Hermitian 

matrix and, consequently, the exponential part of the solution can not be excluded. We have con- 
. jectured the general structure of the solution, both with respect to the quantum number and the 

' order of the expansion. An explicit proof is given for the three leading orders of the asymptotical 

I solution and we sketch a proof for the asymptotical convergence of eigenvectors with respect to 

, norm. From a more practical point of view, we can estimate the required effort for improving the 

' known solution and the accuracy of the eigenvectors. The applied method can be generalised in 

order to accommodate several variables. 

\o : 
> ■ 

00 ' I- INTRODUCTION 

cn : 

I This paper is closely related to one of the famous eigenvalue problems, namely that of a one-dimensional harmonic 
. oscillator. It is common knowledge that if the eigenvectors are required to have continuous second-order derivatives, 
' each eigenvector is expressible as a product of a Hermite polynomial and an exponential term. The corresponding 
, eigenvalues are equidistantly spaced and bounded from below. Another way to state the problem is given by the 
' annihilation and creation operators which directly diagonalise the Hamiltonian. In comparison, the quartic anharmonic 
^ ^, oscillator was solved by Bender and Wu in Ref. A method for finding eigenvalues for anharmonic oscillators was 
I ' created by Meifiner and Steinborn in Ref. |2. A general method for polynomial potentials was introduced recently by 
^ ; Meurice.^ 

• Instead of continuous functions, we consider functions defined only on a discrete, equidistantly-spaced and countable 
I set on R. The obvious advantage of this approach is that it transforms the problem into an eigenvalue problem of 
• • : an infinite-dimensional, tri-diagonal matrix. The corresponding Schrodinger equation is canonically conjugate to the 

. ^ ' Mathieu differential equationi^ Numerical solutions for noninteger orders are naturally obtained by diagonalising the 

^> , very same matrix, see Ref. and the references therein for applications. 

5^ ' In physics, the discretised harmonic oscillator is manifestly realised by the Hamiltonian of an isolated Joscph- 
\ son junction^'^ and the Hamiltonian of the, slightly misleadingly named, superconducting single-electron transistor 
(SSET).^!^ Presently, excited states are seldom considered because of the radical approximations under which the 
Hamiltonian is solved. Even if the excited states are numerically obtained, it is not immediately evident, what 
happens when the coupling is changed. In this article we give an explicit, asymptotical solution for the discretised 
harmonic oscillator which corresponds to strong Josephson coupling in case of the SSET. The same Hamiltonian also 
describes the so-called quantum pendulum, or a particle in a periodic potentiali^^ii 

The corresponding asymptotical eigenvalues have been available for almost fifty years due to the work of Meixner 
and Schafke on Mathieu functions in Ref. (l^ First, by calculating the determinant of the matrix representation 
accurately enough, we continue this expansion by several orders in the coupling parameter. Then we propose an 
ansatz that transforms the matrix equation into an infinite set of algebraic equations and proceed by recursively 
solving these equations. The general properties of the coefficients in the ansatz can be obtained by studying any 
occurring regularities and reinserting these into the solution. Thus, in addition to the eigenvalues, we have successfully 
conjuctured the general form of the asymptotical eigenvectors. In each order of the expansion, the expressions are 
quoted in terms of an arbitrary quantum number, n, whenever possible. The leading terms have been determined and 
rephrased in terms of an arbitrary order, m, too. We find that the eigenvectors are asymptotical solutions of certain 
differential equations, which enables us to obtain further orders in their expansions. 

The only real-valued parameter in the solution is the coupling constant, because all coefficients, both in the eigen- 
values and in the ansatz are rational numbers. As a practical application, the rate of convergence of the solution 
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towards numerically obtained, "exact", solution, can be reliably estimated. In the asymptotical limit, the dependence 
in terms of n and m assumes the form of a simple monomial, at least down to the limits of numerical precision. 

The solutions of order m < 5 are very simple to program and directly apply as numerical solutions of the discretised 
harmonic oscillator. For sufhcicntly small values of the coupling constant the eigenvectors are practically exact and 
thus they facilitate studies which require the structure of the excited states. We have proven, with the help of 
recursion relations of Hermite polynomials, that the first three leading orders of the obtained solution are correct. 
The calculation up to the seventh order should be performed in the future. We also outline an explicit proof concerning 
the normwise convergence of the eigenvectors. The asymptotical nature of the solutions must be stressed. A very 
thorough introduction on the subject has been given been given by Boyd in Ref. [3 

It is justified to ask, is the proposed solution completely new. The answer is, naturally, yes and no. Both discretised 
and discrete harmonic oscillators have been widely studied before. Both cases are related to orthogonal polynomials, 
so the work of Kravchuk^^ and Hahni^ must be mentioned. The discrete harmonic oscillator, where the position 
coordinate is restricted to a finite number of values, is explicitly solved by Kravchuk polynomials as shown by Lorente 
in Ref. Several discretisations of the harmonic oscillator have been previously solved, each giving rise to a 
specific class of generalised Hermite polynomials. Discretisation by an exponential lattice {— S Z}, where 
< 9 < 1, defines the so-called g-deformed harmonic oscillator and generalised q- Hermite polynomials which are 
rigorously discussed by Berg and Ruffing in Ref. 0. For other applications of the q-deformed harmonic oscillators, 
see e.g. Refs. and where other discretisations are reviewed, too. Borzov, in Ref. "2^, considers generalised 
derivation operators as generators of Hermite polynomials and states that the generalised Hermite polynomials either 
satisfy a second order differential operator or there is no differential equation of finite order for these polynomials. 
Many other types of generalisations are also known, see e.g. the multi-dimensional Hermite polynomials of Rosler^i, 
Hermite polynomials orthogonal with respect to the measure j^p exp(— where 7 > — l^SiS^, and parabosonic 
Hermite polynomials^ In the future, it must be established whether the presented class of Hermite polynomials is 
related to the g-Hermite polynomials, if it results from some other discretisation or is it an explicit example of the 
second group of Borzov's categorisation. Complementary results concerning the introduction of distant boundaries for 
the continuous problem are also 'know^M--^ Finally, it should be emphasised that instead of deforming the harmonic 
oscillator, we solve its common-sense discretisation, used especially in numerical calculations. The asymptotical effects 
of the discretisation are explicitly calculated. 

We also briefly consider the abruptly changing nature of the solutions when the coupling constant vanishes. This 
behaviour is evident for both versions of the harmonic oscillator and the Mathieu differential equation. The asymp- 
totical nature of the solutions and the eigenvalues is caused by this divergence. For the Mathieu equation this has 
been well documented, see e.g. Refs. 15 01211271 A more physically motivated approach is given by Bender, Pelster and 
Weissbach in Ref. l28L where e.g. the instanton equation and the Blasius equation are examined. The present methods 
are closely related to these, although we can not carry the calculation as far in the perturbative expansion. This is 
explained by the necessity of obtaining the expansion for the eigenvalues which makes the present problem technically 
more demanding. 

The present method can be generalised in a fairly obvious manner. Other differential equations with analytical 
solutions can be discretised in the same manner if the correct expansions are found for all parts of the solution. An 
easier generalisation is related to multi-dimensional difference equations with harmonic (quadratic) potential terms. 
The existing solution^^ for Hamiltonians of one-dimensional arrays of Josephson junctions become more transparent 
with the help of present formalism. 

The paper is organised as follows. In Sec. m we define the discretised harmonic oscillator and connect it to the 
Mathieu differential equation as well as the continuous case. The solution ansatz and the resulting set of equations 
are reviewed. In Sec. IIIII we quote our conjectures for the general form of the coefficients in the ansatz. We also 
present the explicit values of the leading coefficients. In Sec. II VI we study solving the set of equations which yields the 
asymptotical eigenvectors. Efficient truncations of the set of equations are explained. The effort for improving the 
obtained results with the present method is estimated. In Sec. we prove that the solution satisfies the difference 
equations, at least for the three leading orders. The rate of convergence and the induced asymptotical orthonormality 
are also reviewed. Finally, in Sec. IVII the conclusions are drawn and an outlook of future possibilities is given. 

A final note for those that are only interested in applying these results in numerical and/or theoretical analysis. 
Please review the beginning of Sec. ^] in order to find the correct parameters for the discretised harmonic or Mathieu 
equation. Then proceed to Sec. IIIII and use the given expressions as approximate solutions in Eq. H3U|> . 
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II. THE DISCRETISED HARMONIC OSCILLATOR 



The eigenvalue problem corresponding to the harmonic oscillator is the differential equation for ipix), 

1 (ftp OJ^X^ , , , 

The eigenvectors corresponding to the well-known eigenvalues, 

A„=a;(n+l/2), (2) 

where n = 0, 1, 2, . . ., are given by 

M^)=AnH„iOe~^'^\ (3) 



Here ^ = y/uj x, An is a normalisation factor, and _ff„ is the Hermite polynomial of order n. The Hermite polynomials 
are solutions of the Hermite differential equation 

y" - 2xy' + 2ny = 0, (4) 

where n = 0, 1, 2, . . .. For our convenience, we write the polynomials, given by Rodrigues' formula, as 

fe' 



^ fc=0 



where k' :— \n/2\, i.e. k' = n/2 if n is even and k' — (n ~ l)/2 if n is odd. The quantity k' proves to be extremely 
useful in further analysis. The Hermite polynomials satisfy the recursion relation 

Hn+iiO = 2eH„(0 - 2nff„_i(e). (6) 

Many of the generalisations of the Hermite polynomials boil down to a generalisation of this recursion relation>iSiiLi2i24 

The discretised version of Eq. |^ is obtained by restricting the values of x onto an evenly spaced, countable subset of 
R. This corresponds e.g. to the discretisation of charge in case of a Josephson junction or a SSET. Only the constant 
nearest- neighbour coupling is retained which yields a tri-diagonal matrix H{xq) with non-zero matrix elements 

Hjj{xo) = \uj'^{j - a;o)^ i?j+ij(.To) = Hjj+i{xo) = ~\. (7) 

Here the parameter xq G I^^' 5] ^^"^ displacement of the origin with respect to the matrix element j = 0. All 
eigenvalues of H{xo) have been translated by —1 in order to simplify the diagonal matrix elements. The standard way 
to write the Hamiltonian of an inhomogeneous SSET is obtained from Eqs. (7.36) and (7.39) of Ref. ^3 and rephrasing 
it in terms of the number operator for Cooper pairs yields the matrix 

H^f''^\No) = Ecij - Nof, Hf^fPiNo) = {xo) = ^^E,{e), (8) 

where Nq is the number of Cooper pairs which minimises the charging energy, Eq = (2e)^/2CE is the unit of charging 
energy, and Ej{9) is the effective Josephson energy which depends on the total phase 6 across the SSET. Consequently, 
we solve the Hamiltonian of SSET if we find the eigenenergies and eigenvector for the discretised harmonic oscillator 

withw = (2£:c/£^J(6'))^''^• 
In the following, we are searching for eigenvectors with finite Euclidean norm, i.e. 



00 



E 

j=~oo 



|^,f<oo. (9) 



The existence and uniqueness of such solutions follows from the generalisation of the Gershgorin eigenvalue theory 
by Shivakumar, Rudraiah and Williams in Ref. Is^- First the number of eigenvalues of H{xo) on a given interval 
can be shown to coincide with number of eigenvalues for a finite-dimensional truncation of the matrix, _ff*^^^(a:o), if 
the dimension N is sufficiently large. A sufficient condition for this is that the difference between ordered diagonal 
elements exceeds 2 x | — i | — 1. When x ~ 0.5, this property is obtained more easily for even values of N. Furthermore, 

they prove that, for finite values of n, the eigenvector ipn^'' of H'^^^{xq) tends to the corresponding eigenvector of 
H{xq) when N — s- 00. 
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We now establish the connection between the H{xo) and the Mathieu differential equation^ 

fl + {a-2qcoai2v))y^0, (10) 

where a is the eigenvalue, also known as the characteristic value when the solution y has period of tt or 2tt. We follow 
the derivation of Shirts in Ref. and use Floquet's theorem to obtain 



y ~ exp(ih'v)P{v) = exp(ww) C2k exp{2ikv) (11) 



where the Fourier expansion of P{v) has been inserted. This corresponds to the matrix equation for the coefficients 
C2fc compactly written as 

C2k-2 - V2kC2k + C2k+2 = 0, (12) 

where V2k — [a — {v + 2k)'^]/q. This is identical to the discretised harmonic oscillator Hamiltonian H{xq) with an 
eigenvalue A after identifications 

1^ = 2x0, k^j, g = 4/w^ a = 8A/w2. (13) 

Thus all results obtained for the discretised harmonic oscillator also hold for Mathieu functions Hll|) with parameters 
given in Eq. \i'6l . For xo = and Xq = ±i the solutions of H{xq) can be chosen to be even or odd with respect to j. 
This corresponds to writing P{v) in terms of sines and cosines. Special attention must be given to the even solutions 
of H{xq = 0), where the resulting equations in the matrix representation read 

-^Pi/V2 = A2„Vo, (14) 
-W^ + wVi/2-V'2/2 = A2„i^i, (15) 
-V^,-i/2 + c<j2j2^j./2- V,+i/2 = \2n4>3- i>2. (16) 

The eigenvalues for xq = correspond to characteristic values {a2n{q),b2n{q)}i while the case xq = ±i is linked to 
{a2n+i((/), 62n+i('?)} as defined in Ref. 

The asymptotical expansion of the eigenvalues corresponding to the limit q — > oo or w — > was obtained by 
Meixner and Schafke in Ref. The derivation of the eigenvalues is based on the three-term recurrence relations for 
the Mathieu functions and the requirement that the norm of the error in the eigenvalue equation vanishes faster than 
a specific power of uj. Meixner and Schafke quote the asymptotical characteristic values of the Mathieu equation up 
to and including the order cj^ in Theorem 7 in Sec. 2.3. Some error estimates for asymptotical expansions of Mathieu 
functions by M. Kurz are given in Ref. 27, Because the Mathieu equation is also the Schrodinger equation of the 
quantum pendulum or a particle in a periodic potential, it has been studied independently in physics, too^S^ii^ 
Especially, the same general expansion for eigenvalues and several further terms for the ground state energy were 
obtained by Stone and Reeve in Ref. 

In this limit, we can write the eigenvalues of H{xq) as 

oo 

where w — > 0, and 

m' 

Ai"^ =5]A|:;;)fi™+2(^-™') (18) 

fe=0 

with n :—2n+l and m' = \ rn/2\. This structure is identical to that of the Hermite polynomials Q, if one identifies 
n with By Ref. 0, the eigenvalues (^0)1 do depend on xq, but this dependence decreases exponentially as ^ 0. 
The maximal difference is given bjii^ 

A„(xo = ±i) - A„(xo - 0) ^ (-l)"So(l - exp(-8/c^), (19) 



where Bq and Bi depend on n but not on uj. 
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This allows us to write the eigenvalues of H{xq) 



as 

A 







Lo'^di Lo^dr-, ui^dg 


up dj 


u^ds 






211 


217 223 227 

Lu^^di3 oj^'^du 


233 

w^^dis 


240 

w^'die 


247 


257 


261 


269 272 


279 


287 



251 

where the coefhcients d^ read 

d2 = ri^ + 1, 

ds = ri"^ + 3n, 

d4 = Sn"* + 34ri2 + 9, 

ds = 33n^ + 41071^ + 405n, 

dg = 63n^ + 1260n* + 2943ri2 + 486, 

dy = 527n^ + 15617n5 + 69001n3 + 41607n, 

ds = 9387n^ + 388780n^ + 2845898n'' + 40218847i2 + 506979, 

dg = 175045n^ + 9702612n^ + 107798166n5 + 288161796fi^ + 130610637n, 

dio = 422565^1° + 30315780n^ + 480439190n'^ + 2135766820n'' + 2249346285^^ + 238353840, 

dii = 4194753ri" + 379291385n^ + 8186829426n^ + 55529955498^^ + 110241863469n^H 
41540033277n, 

, n 

dl2 



(20) 



^13 

di4 = 

di5 = 



41540033277n. 

^ 10645960ni^ + 1187264199^^° + 33678377895n** + 327725946398/1^+ 

1081358909790n^ + 94007705503571^ + 88258370067 
= 440374207^13 + 59495737574n" + 2155821044201n^ + 28738150160500n^+ 

144821249264769^5 + 236410740537606^^ + 78243613727607fi 
= 578183175^1'* + 93209584104ni2 + 4215683624295ni" + 74269604367684n^+ 

537905750769429n^ + 145676730601375271" + 11057115504106537j2 + 94839535889532 
= 12308013927ni''^ + 2337227706555771^ + 12943725324367577^ + 292850645568409577^+ 
29119560960614085n^ + 120372998803922241n-'^ + 17092192064940274577^+ 
5131634402399008577 

= 53003912615977*^ + 117243302735480771* + 782309396142565277*2 + 22204381081902685677* 
292495292113002519477* + 1738031526802826522477^ + 4085166941152660098077*+ 
2798355147033036578477^ + 2235152520630714879 



We obtain the terms for orders % < m < 11 by exploiting Eq. (|18() when explicitly evaluating the determinant of 
Eq. {Tj). As a first step, setting = halves the dimension of the tridiagonal matrix. Next, by translating one of 
the eigenvalues close to zero by substracting the known expansion of this eigenvalue, the determinant becomes an 
essentially linear function of the chosen, translated eigenvalue. The next unknown term is inserted as a parameter and 
the determinant is calculated for several values of w, preferably in the form {2^'^a;o}|=o ^ ■ This choice lets us separate 
the leading correction and the subsequent corrections. In order to obtain the terms d2fc and d2k+i, we must correctly 
determine all eigenvalues A„ when n < k. SufRcient accuracy is guaranteed by using the high-precision numerics of 
Mathematica software. The method for obtaining the orders 777 > 11 requires explicit knowledge on the properties 
of the eigenvectors and the discussion is postponed until the end of Sec. 11111 

The asymptotical nature of the expansion means that for each value of uj and 77, there exists and optimal order 777 
which minimises the error in the eigenvalue, i.e. the function 



AA(w, n, m) 



m'=0 



(21) 



with respect to 777. The exact eigenvalue A„ exists and is finite for all non-zero values of uj according to the Sturmian 
theory of second-order linear differential equations, see e.g. Ref. |5j. In order words, for sufficiently small values of w, 77 
and 777, the error is dominated by the first omitted term, i.e. AA(a-', 77, 777) ~ jAi™"''^'' jw^+i. Because the asymptotical 
eigenvalue is divergent, it surely crosses the exact eigenvalue when uj is increased, but this occurs outside the range 
of asymptotical convergence. Similar asymptotical convergence should be observed for the asymptotical eigenvectors, 
too. Assuming ipl^''^"^ corresponds to the asymptotical expansion of the eigenvalues up to and including order w™, 
we expect error in the norm to behave as 

||^(m..„)_^(.„)|j^C'(77,r77)c.™, (22) 

where a; — s- and C(77, 777) is a simple function of 77 and 777. Although this has not been proven, Eq. H22|l appears to be 
correct and we will ultimately give an approximate expression for (7(77,777), too. Outside the regime of asymptotical 
convergence the error (|22|l approaches V2 as the asymptotical solution becomes orthogonal to the exact one. 
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Next we show that the discrete eigenvalue problem Eq. Q is a meaningful asymptotical limit of the continuous 
harmonic oscillator equation The problems are identical in the leading infinitesimal order when co is infinitesimal, 
but the limit a; ^ is subtle. As long as a; > 0, both the eigenvalues and eigenvectors of the discretised problem tend 
to those of the continuous harmonic oscillator with this lu. For lu — the continuous problem becomes abruptly the 
free particle Hamiltonian with solutions 

i^^=o{x) = e"'^ , A^=o = fcV2, (23) 

where k is the standard name for the wave number. Simultaneously the discretised problem becomes the well-known 
nearest-neighbour chain with eigenvectors and eigenvalues 

^^ = {e^fcO-.o)}^.^ A„=o = -cos(A;). (24) 

For sufficiently small values of k we have A^^^o ~ —l + k'^/2, in agreement with Eq. (|23|l . In contrast, we are interested 
in the bound-state solutions of Eq. 0} and those eigenvectors of the discretised problem that can be uniquely related 
to these continuous solutions for u; > 0. 

The harmonic oscillator is discretised by restricting the values of x onto a countable and evenly spaced subset of 
K. The lowest-order central approximation for a second-order derivative is simply 

Assuming ipi^) to be real-analytic allows us to write the numerator of the right-hand-side as a Taylor series 

°° r,i,2k J2kj,(^\ 

^{x + h)- 2i.ix) +.pix-h)=Y: 1^^^- (26) 

If h is infinitesimal and as the derivatives of are finite in all orders, the only remaining term is h?%l)"{x). Thus, in 
the lowest infinitesimal order the discretised eigenvalue problem gives a second order differential equation 

1 cPipx w^a;^ 



2 dx^ 2h^ 



= /i-^(-l + A)V'. (27) 



which is identical to Eq. JQ) apart from the constant —h^^ and the redefinitions lu ^-^ ui/h and A i— > Xjh?. The 
discreteness of the problem can also be varied by rescaling the value of lo. Thus, instead of decreasing the size h of 
the steps, we set h = 1 and let ui ^ 0. From Eq. H27(l we see that asymptotically the eigenvalues and eigenvectors 
have the form A„ ^ — 1 + uj{n + 1/2) and V's ^ as expected. 

We have already pointed out that the matrix H{xo) in Eq. iQ can be derived from the Mathieu equation. The 
underlying reason for this is that the problems are canonically conjugate. Inserting the full expansion Eq. (|26|l into 
Eq. (|27|l yields an obvious differential equation in with respect to x. The canonical transformation id/dj v and 
j — > —id/dv preserves the eigenvalues and produces the differential equation 



\k~2k 



2 dv^ (2fc) 



= X^i- (28) 



Noticing that the sum is equal to cos(5) and setting v :— {ii + n)/2, we obtain the canonical form of the Mathieu 
equation with parameters given in Eq. 

After these important preliminaries, we are able to proceed towards the actual solution for the discretised harmonic 
oscillator. In order to treat eigenvectors of all matrices H{xo) on an equal footing, we replace the index j by x := J — xq. 

(n) 

For arbitrary values of xq and j the new index x becomes a continuous one on R. We thus obtain functions tpx , 

(n) 

where n is the state index. We propose that these functions ipx are real-analytic and that they give the eigenvectors 
of H{xo) asymptotically, i.e. 

- {^j-loi.^-oo (29) 

when u ^ 0. The problem tends to the continuous one in the lowest (infinitesimal) approximation in lu. Thus it 

is reasonable to assume that the lowest-order approximation for the solution functions is given by ip^'' ^ ^n{x) as 
w ^ 0. 
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The general form of the asymptotical solution of the discretised harmonic oscillator now reads 

Coo CXD \ k' oo 

fe=l l=k / k=0 1=1 

where a^"'' and /J^"' are constants to be determined. The solution to the continuous case yields a^^^l = — 1/2 and 

(n) (n) 

/3j, { = 1. We are free to normalise the solution so we can choose /3q ^ = for ^ > 1. 

(n) 

The main point of introducing the functions Tpx that they transform the difference-equation-type eigenvalue 
problem corresponding to the discretised harmonic oscillator into an infinite set of algebraic equations for each value 
of n. The eigenvalues (|2U|) appear as parameters and they are required in order to solve the equations for the sets 

of coefficients {a^."'*} and Fortunately, the equations uniquely determine every single coefficient. Because the 

expansion of the eigenvalues is asymptotical, the meaning of the full solution to these equations must be determined 
later. 

In practice, we need a suitable truncation of Eq. H3U|) and thus we define an (unnormalised) approximate eigenvector 



(31) 



where ipi^'"^^ contains only those terms with I < m. The definition of ijjn'^°^ obviously coincides with the continuous 
solution at xq. In numerical calculations, and always for even values of m, we must truncate the eigenvector with 
respect to j, by setting {^p'^^^''^°'')j = for components |j| > jo with a sufficiently large value of jo- 

We now give the infinite set of algebraic equations corresponding to the transformation of the difference equation 

(n) 

when the solution functions i/ji are substituted into the eigenvalue equation. Rearranging the terms, we find that 
each equation can be written in the form 



i - ,/,(») 



V'r(-A„+c.V/2), 



(32) 



where x — j — xq. Inserting the general ansatz l|3U|) into Eq. (|32|l expresses the equation in terms of ^ and lu. The 
exponential part of the ansatz on the right-hand-side canceled simply by substracting the corresponding exponent 
from those on the left-hand-side. This yields an equation 



oo oo 



-p ( E E 4?^-'=+'-M(- - 1)^^ - E E [v^i- ~ i)r'^''''^) + 

\k=l l=k / fe=0 1=1 

/ OO oo \ k, OQ 



exp EE4';^-'^'-M(-+i) 

l=k 

k 



2k ^2fci 



k=0 1=1 



' \i{7n) 
m=0 



(33) 



k=Q 1=1 



These equations are then expanded as functions of x and uj as the resulting equations are easier to solve. The equations 
must hold for all values of the linearly independent variables x and uo so each equation must be solved separately. For 
the purposes of generality, it would be preferable to expand with respect to ^, but the resulting equations are much 
more difficult, both to obtain and to solve. Nevertheless, the obtained solution can be inserted into to these equations 
in order to show that the results are correct. This will be done in Sec. IVI 



In order to obtain the eigenvector i/; 



we must solve and satisfy all equations corresponding to 



n'\'2'rn' — 21' \m' + k' 



(34) 



This is of course done recursively, by inserting the known part of the solution and solving for the next level. In order 
to connect Eq. (|34|l with the order of the solution, we state that the equations corresponding to a fixed value of m' 
uniquely determines the coefficients with I — m' . 

After a while, one starts to see regularities in the coefficients and attempts to express these in a functional form. 
We have been able to find rather general expressions for the coefficients. This means that the coefficients have been 
expressed in terms of n and the order of the expansion, whenever this is possible. We have conjectured the general 
form of the terms which means that we know how far away we are from obtaining further terms. 
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We have found that the functions asymptotical solutions of the differential equation 




(n,m) 

fe=0 



^AWc.'^^i"-™) (35) 



in a specific sense. After all derivatives have been carried out, the terms multiplying the common exponential factor 
cancel up to and including the order w™. If the solution i/'i"'™"^-' is known, we obtain an explicit differential equation 
for the exponential part w'"~i/m(^), the correction to the Hermite polynomial uj^^^^gmiO ^^'^ the energy eigenvalue 
Xt\ In case of the ground state and the first excited state {n — 1), the condition gm{Cl = renders the problem 
solvable. For n > 2 we must insert the ansatz (|3UII in order to obtain the solution. 

The results are, naturally, in complete agreement with those obtained by using the difference equation. We are 
using just another representation of the original problem. Equation l|35|) enables us to obtain the solutions for fixed 
values of n up to relatively high orders with respect to powers of lo. Thus we can both extend the general expression 
for the eigenenergies in Eq. I|20() and those for the coefficients in the exponential part of the solutions. For the ground 
state energy, we find the terms beyond order oj^^ to be 



363372562420411197^1^ 6258692522467212813 227867608383920243815 w 
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279 2^3 2^8 

4372199488222446620121 w^o 352807992522448740907163 w^i 7465886451386334274097895 



292 2^8 21°^ 

330752735437897260202410959 cj^s 7654237307570898665851927581 



2107 2111 
1477812451863756884805687589129 w^s 37132718819258763418452357390369 w^e 



2118 2122 

1939848955425261040700592191917783 cj^t 52598573101029275526869814635336865 w^s 



2128 2131 

5914101566562517015636997146651378649 172129355454985486683952198830698506149 w^o 



2137 2141 
10362392343003738344189045786484697182753^31 

^TiS +0{uj^^). (36) 

The corresponding asymptotical eigenvector contains 31(l + 31)/2 = 496 linearly independent terms. The coefficient 
of uj^^x^ in the exponential part reads 

5207328980459439428858189871778019425519567564728193 

2765292404617797269550429065808396826741571584 ' ^ 

The general solution is given in the next section. 



III. THE GENERAL SOLUTION OF THE DISCRETISED HARMONIC OSCILLATOR 

For reasons of completeness and easy accessibility some of the definitions will be repeated in this section. The 
j^th Q].(jgj. solution function corresponds all terms up to and including Z = m in Eq. (|30|) . The asymptotical 

expansion of the eigenvalues A„ is given in Eq. (|20|l . The state index n determines two expansion parameters, 

fi := 2n + 1 and k' [n/2j, (38) 

where k' = n/2 if n is even and k' = [n— l)/2 if n is odd. The gamma function T{x) is the generalised factorial with 
the defining property xT[x) = T{x + 1). We need the values for integer and half-integer values which read 

r(A:) = (fc- 1)!, ^(fc + i) = 2-^^/¥(2A:-l)!!, (39) 

where the double factorial fc!! for integer values of k is given by k{k — 2) x • • • x (1 or 2). The coefficients in the 
Hermite polynomials simplify to 



^(„) _ (_i)fc'+fc22fc+(i-(-i)")/2„! 

(2/c+ (1 - (-l)«)/2)!(fc' - /c)!' 
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A convenient normalisation for the eigenvectors is obtained by requiring that 

^H^^(i-(-TO^ (41) 

i.e. ~ 1 for even values of n and ~ ^ for odd values of n. Also bear in mind that 



ai':/ = -l/2, and /3^7(>i) = 0. (42) 



= -1/2, /?H = 1, and (3^-\ 
Under these constraints we have conjectured that the general form of the coefficients. In the exponential part, 

/ OO CO \ 

exp EE4:^-"^^^M' (43) 



the coefficients can be written as 



a 



\k=l l=k 



k,k+l 



E4'i+z"'- (44) 



;'=o 



In) 

Please note that if the coefficient a\i are written as polynomials in n instead of n, the signs of the corresponding 
expansion coefficients 6^1 appear to be given by (—1)*^. An efficient way to write these coefficients is given by 

where Q{k, l,2k ~ 5 + I') and Q{k, I, k — 2 + I') are polynomials in k of orders 2k — 5 + I' and k — 2 + I', respectively. 
An important consequence of Eq. H45|l is that regardless of the values of / and n we have 

,1^^ = -1/4. (46) 

The explicit expressions for the seven leading coefficients have been obtained and they read: 



a 



(„) _ (_i)fc2^-^^r(fc + i/2) 

fc(2fc- l)2r(A:)0F 



kk 



^{n) _ / -,xfe^-2-2fc M , r(fc + 1/2) 



'=■'=+1 ^ ' \k T{k + k J ' 

- (-l)^-2-^-^^^(n+^I|^^[(3 + 52fc + 40fc2) + (9 + 12fc)n2]) 
4") = (,i)fc2-9-2fc/(_i + 7fc + 5fc2) + (3 + 4fc)n2+ ^^(^ + ^ 



24r(fc + 3)V^ 

[(243 + 1119fc + 1928A:2 + 1376fc^ + 320fc'^)n + (33 + lOlfc + 104fc2 + 32k^)n^] } , 

4"' - (_i)fe2-"-2fe 1(53 ^ ^20fc + 136^2 + 40fc3)n + (37 + 72fc + 32^2)^3/3 + ^(^ + 1/2) 
'''''+^ [ 48r(A: + 4)0r 

[(-2612925 - 5292132yfc + 10675063^^ + 36766856/c^ + 40148416A:'' + 21300608fc^ + 5544448/c^+ 
565760fc^)/315 + (11070 + 60044fc + 130810fc2 + 142112fc3 + 81280fc* + 23168fc^ + 2560fc^)n2 + 
(585 + 2288fc + 3585fc2 + 2696A:^ + 960fc'' + 128fc^)n'^] } , 
4"^_^5 = (_i)fe2-2o-2fe |(_5i87 _ 672fc + 6580^^ + 7684fc3 + 3164fc'^ + 452fc^)/3 + (1214 + 3744fc+ 

3^^4/o , r(fc + l/2) 



4080fc^ + 1968fc'^ + 320fc'*)n'' + (345 + 808fc + 576fc^ + 128fc-')n73 



48r(A: + 5)V^ 

[(740893230 + 3944788389A: + 9627147810/c^ + 14943869467A:^ + 15287941200fc'*+ 
10116675072fc^ + 4238798592fcS + 1079918592A;^ + 152076288A:^ + 9052160fc^)n/315+ 
(1825740 + 11037114fc + 27955236fc2 + 37919062fc3 + 30169312fc'^ + 14491648fc^ + 4122880A:^+ 
636928fc^ + 40960fc*^)nV3 + (85050 + 381087fc + 729798fc2 + 752369fc3 + 447024^;-*+ 
152576fc5 + 27648fc^ + 2048fc^)nV5] } , 
4"^^g = (_j^)fe2-26-2fc 1(378033 + 496368fc + 786528fc2 + 710816fc3 + 339904fc'* + 79552/?^ + 7232fc^)n/3- 
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(69714 + 241312fc + 303392fc2 + 177696fc3 + 49408fc^ + 5120k'^)h^ /3+ 
(17217 + 45360fc + 40960A:2 + 15360/0^ + 2048/c'*)nVl5 + ^''^ ^ 



18Or(A: + 6)0F 

[(-24640192386810- 105728184475128fc- 155775948330744fc2- 74654535511116fc3+ 
74660144680858fc^ + 156803802177352^^ + 134434233033760fc^ + 70722102090816fc^+ 
24590691451392fc^ + 5680345583616fc^ + 839668527104fci° + 71921254400fc"+ 
2714009600fc^^) /9009 + (22093103970 + 162201234402fc + 504160865145fcV 
882850470198fc^ + 986932878421^" + 745434338828fc^ + 388089936864^^^+ 
138972684672fc^ + 33504543744fc** + 5179637760fc^ + 462565376fc^° + 18104320fc")nV21+ 
(152041050 + 991922940fc + 2784482730^^ + 4353707520fc3 + 4203836660fc"+ 
2632731680fc^ + 1088777440^;*^ + 294912320fc^ + 50245120fc* + 4874240fc^ + 204800fc^°)n'*+ 
(2606310 + 12799746fc + 27798345fc2 + 34245070^^ + 26181505^;^^+ 
12857468fc^ + 4055200^^ + 792320fc^ + 87040fc* + 4096fc^)n^] } . 

Thus, for an arbitrary order m, we can obtain the expressions for coefheients corresponding to {^™-i^2m-i }6^^ 
Furthermore we find 

a^"g^ = -(505549159n+ 177209155^3 + 8289645^5 + 40329n^)/237- 

(-2741702 + 1224882571^ + 1518052n* + 26073n^)/23^ 
a["^ = -(-840819020949+ 1419128841068^2 + 221074444682n'' + 6195597884n^ + 21259875fi^)/(3 x 2^'^)- 

(131878584971 + 45 9 3 8 9 2 5571^ + 2971811l7i^ + 335617fi^)/23^ 
4"g' = (13125727618771 + 3784309918771^ + 13230464977i5 + 4456305n^)/2"''+ 
(48228434 + 9395984577^ + 87877007i" + 110661riS)/23'*. 

The exponential part H43|l is now completely determined up to the ninth order, i.e. known for arbitrary values of n 
for terms with ^ < 9. 

If we exclude the dependence on 2^^^^ and also that given by the gamma functions in the coefficients, we observe a 
very distinct regularity. The dependence of the leading power k in each polynomial sequence starting from nP for a 

(n) 

fixed value I in ^_|_; and going upwards by one both for I and the power of n is so far always given by 

{5(u')h' = {%o)/(4'' (/')!)}■ (47) 
The initial values in cases T < 5 are given by 

{a(/,o)}f=o {!> 1/4, 5/48, 5/512, 221/96768, 113/786432}. (48) 

This dependence is by no means proven but it corroborates our choice for the prefactors in Eq. H45|l . 

The coefficients {/?^"^} determine a set of new polynomials, where the coefficients multiplying the powers of ^ 
depend on ut. In the limit uj ^ these polynomials tend to the Hermite polynomials. They are, unquestionably, a 
new class of generalised Hermite polynomials. They are defined as parts of the eigenvectors of a Hermitian matrix. 
Because the exponential part is rather complicated, the measure, with respect to which they become asymptotically 
orthogonal, is necessarily a complicated one. It depends on both of the eigenvectors, i.e., it is not a measure in the 
classical sense at all. No simple recursion relation for the polynomials is yet known, and we do not know, whether they 
satisfy any differential equation of finite order. This means that they could be an example of the second category of 
generalised Hermite polynomials as defined by Borzov in Ref. Such discussion is beyond the scope of the present 
study and we will concentrate on simpler properties of the polynomials. 

Our generalised Hermite polynomials are defined as 

^^n(6:=E4"^r+^''=-'='\ (49) 

fe=0 

where the modified coefficients are given by 

oo 

1=1 

Because the generalised Hermite polynomials i?" (^) fix the nodes (zeroes) of the functions t/iI"-* , it is equally important 
to obtain correct polynomials as it is to obtain the correct exponential factors. 
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We conjecture that the general form of the coefficients {P^i^} reads 

= E ( * E (pS + [(1 - (-1)")/2]P-S)fc' 1 {k'f, (51) 

l'=0 \ 1 = 1 

where and arc constants. Additionally, p^l^- = when I = 2{l — 1) — V or I' = I — 1. This expansion with 
respect to k and k' shows that even and odd values of n should be treated separately. 

Some general properties of the coefficients /o^'^- and Pj'^- have been gleaned. The recurring appearance of the factor 
(lOfc' — k) is by far the most striking of the observed regularities. This factor may, in time, explain some properties 
generalised Hermite polynomials. We conjecture that 

E (/^I'ic^-D-n.o-r^'^'"'^""""'"') = k^-'-'^^'ilOk' - fc)'-^-^"<'P(2mo,0, (52) 

l'=0 

where 2mo < / and P(2mo, denotes a (2mo)*'^ order polynomial in k and k'. Similarly, the diff'erence between even 
and odd values of n corresponds to 

1-2 

E (p|';2(-3-n^o-r^''"'""°"''"0 =^'"'""°^^°^'~^^'"'"'"°^(^"'°'^^' (^^^ 
l'=0 

where 2mo < / — 1 and P{2mo,l) again denotes a (2mo)*'' order polynomial in k and k'. 

In the leading and ncxt-to-lcading orders the polynomials P{0,1), P{2,1), P{0,1), and P{2,1) have been explicitly 
evaluated. Thus, we define the quantities 

B{l,2l-2).- 48,_i(^_^^, , (54) 

and 

~ 5x48'-i(/- 2)! ~ 2)658(fc')' + (402 - 126/)^'^ + (8/ - 31)k% (55) 

which has been confirmed up to the sixth order, i.e. I = 6. Similarly, the leading differences give rise to the quantities 



B{l,2l-3) 

and 

B{l,2l - 4) := 5^4gi-i(^_3)| [(2632; - 2576)(fc')' + (1470 - 504/)fc'A;- + (32/ - 145)^;^]. (57) 

The leading terms arc very similar, but also the ncxt-to-lcading terms B{l,2l — 3) and B{l,2l — 4) share several 
common features. Most importantly, the /-dependence in the polynomial section is identical, apart from a factor of 4. 

Below, we give the explicit values of the coefficients /3j."^ in cases 2 < / < 7. It is convenient to separate the even 
and odd values of n, because the correct expansion parameter appears to be k' . Please note that these expressions 
automatically yield /3q"^ = for / > 1. 

0fV = [(3fc-fc2) + (10fc)fc']/48, 

f3^^P = [(855fc - 64fc^ - 14A;^ + 5fc^) -t- (784A; + 48A;^ - lOOfc^)^' + (1316/c + 500k'^){k'f] /23040, 

/jg'+i) = p'^^p + [{249k + 49k'^ -20k^) + (532k + 200k'^)k'] / 11520, 

= 5(4, 6) + B(4, 5) + [(371385fc - 203498/0^ - 12129fc-'' + 1438fc'^) + 

(1110698fc + 102042fc2 - 26252k^)k' + (496932A; + 93984fc2)(fc')^ + 560200A;(A;')^] /23224320, 



4fc'-i(10fc'-A:) 
48'-2(/-2)! 



1-2 



(56) 
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/3^2^^'+i) ^ + Bii, 5) + 5(4, 4) + [(67680fc + 12347fc2 - 2602fc3) + 

(108544A: + 33762fc2)fc' + 114456A:(fc')^] /3870720. 
(3\^P = 5(5, 8) + B(5, 7) + [(278751375fc - 202014918fc2 + 35222268fc3 + 4026748^"^ + 28158fc5 - 9944^^) + 
(713250468fc - 281790420fc2 - 61452368fc^ ~ 196176fc'* + 209856fc^)fc'+ 
(1105743252fc + 198178852fc2 - 10630680fc3 - 275344fc'')(fc')^ + 

(319197168fc + 81282336fc2 - 22799744fc^) {k'f+ (271672512/c + 148408976fc2) [k'f] /22295347200, 
/3^^5^'+i) = f3l^P + B(5, 7) + B{5, 6) + [(119817225fc - 23468037^^ - 12060122^^ - 330312^-* + 100198fc^)+ 
(436319556fc + 103769756fc2 - 5886336^^ - 1921112fc'')fc'+ 

(321608148fc+ 118383396fc2 - 904080k^){k'f + (224147856fc + 120486304fc2)(fc')3] /11147673600 
(3\^P = 5(6, 10) + B(6, 9) + [(134035780725fc - 166751340588fc^ + 39327194883^^ - 2269605874^"*- 

477614210fc^ - 19226552fc6 + 221782fc^ + 484fc*) + (413990823078fc - 217584747090fcV 
+22678956764yfc=^ + 8841166604^"* + 479019924fc^ - 6549884fc^ + 566984fc^)fc'+ 

(526339688532fc-155591533528fc2-55003198072fc3-3518713436fc'' + 85514000fc^-24903296fc^)(fc')^ + 
(556945898088fc4- 131085561976fc2 + 2790556248^^ - 280060176fc"* + 424040144fc5)(fc')^+ 
(116760015552fc + 34523271136fc2 - 5865150192fc3 - 3338174576fc"*)(fc')''+ 
(79966766400fc + 46102886720^^ + 10162787360fc3)(/c')^] /11771943321600, 
= ^(2fc') ^ ^^^^ ^ ^^^^ ^ [(34460588160fc - 26910050283^^ + 72069996fc3 + 1282383895fc"*+ 
103465570fc^ - 948002fc6 - 316976fc^) + (216801198648fc- 21671791146fc2- 
18471533106fc^ - 1699322576fc'' + 78651584fc^ + 5865684fc^)fc'+ 

(345295895928fc + 96181762100fc^ + 1478206984fc3 - 1409258180fc"* + 60822872fc5)(fc')^ + 
(160052617776fc + 64547633160^^ + 3204992824fc3 - 1898232512fc'')(fc')^+ 
(83156900448fc + 47130830560^^ + 10276562912fc3)(/c')4] /5885971660800, 

p'i^P = 5(7, 12) + 5(7, 11) + [(21167446950775125fc- 34318046368345140fc2 + 13674300462898392fc3- 

1352901404372446fc'' - 2843855572731fc^ + 11311875159790^^ + 704407032828 
12949326156fc* - 177366189fc^ + 37677640fc^°) + (59570630372492640fc- 
60644270495554704fc2 + 10066261 151648252^^ + 104602336760652fc'' - 246415137367020^^- 
20207362771548fc'^ - 460168946016fc^ + 2604105504fc^ - 2483040560fc^)fc'+ 
(99669485611466412fc- 39020273844707836fc2 + 2200698814542984^^ + 2070713072954600fc''+ 
212428368788100fc^ + 5622413614220fc'^ + 82814211480fc^ + 70784553840fc®)(fc')^ + 
(82488078028378080fc- 18442822328400480^^ - 9100818756007520^^ - 964844434165920/c'^- 
24156442527360fc^ - 3013682511840fc^ - 1116228072960fc^)(fc')^ + (66588038149135200fc+ 
18345507366303440fc2 + 1191268975557840fc3 + 10518809509520fc'' + 42474114642960fc'^+ 
10244315921840fc^)(fc')'' + (10839030004200960fc + 3516288982521792^^ - 363895953410496^^- 
304221200739456^"- 51610667908800fc5)(fc')^ + (6218212960526208fc + 3705496740373376fcV 
898601964676416fc^ + 110684037464000fc"*)(fc')^] /1542595452862464000, 
^ ^(2fc') ^^^^ ^ -^Q^ ^ [(-460686821541975fc - 1941941074537755^^ + 366877584331212fc3+ 

41494582964306/c'' - 9965970910165fc^ - 1112954021925fc^ - 34972438722fc^ + 958101144fc*+ 
3798795fc^) + (9145976126266080fc- 3823250702059872fc2 _ 191143081971676fc3+ 
173415983605096fc"* + 24948282593576fc^ + 815890000796fc'^ - 41300603344fc^ + 1974092120fc*^)fc'+ 
(19045045703842332fc- 607795420829248fc2- 1422456568355676fc^- 195873935369616fc"*- 
2655255816252fc^ + 611948653316/c^ - 98001423520fc^)(fc')^ + (18941405236672032fc+ 
5863460843089248A:^ + 345511721120640fc^ - 55753844615456fc* - 654723253184A:^+ 
1792716525600A:^)(fc')^ + (6296099301099168A: + 2692631923242160fc2 + 232744611197904fc3- 
59049744898144fc'' - 14606140268720fc5)(A:')'' + (2605202959125888fc + 1525593370591680fcV 
365590616623232fc3 + 44670372947200/0") (fc')^] /257099242143744000. 

In combination with the exponential parts these coefficients determine exphcit, analytical expressions for solution 
function ^i"'^' for arbitrary values of 

It must be re-emphasised that Eq. I|30|l is an asymptotical solution. Two partially overlapping reasons for this 
behaviour must be stated. First, the solution depends on two independent length scales, i.e. x and cj, and second, 
the coordinate transformation x i-^ — ^/ujx is singular at a; = 0. These points are rather extensively covered in 
Ref. The eigenvalues are asymptotically exact for even values of n at = and, probably, for odd values of 
n at ccq = ±i. Because the error decays exponentially in 1/w, this dependence on xq vanishes much before the 

asymptotical behaviour of Eq. if^ . i.e || V'i™''^'''' — V'"^"''!! ^ C(n,m)a;™, appears. 
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Comparison against numerically obtained eigenstates allows us to give an approximate expression for the function 
C(n,m). The validity of the calculations is limited by the numerical precision, i.e. to norms of the order of 10~^^- 
10^^^. We have employed the reliable diagonalisation routines of Matlab software for this purpose. We have studied 
eigenvectors up to n « 40-50 and the corresponding asymptotical solutions '^i'"'^"'' up to the fifth order. A reasonable, 
order-of-magnitude estimate for the error in the Euclidean norm, when n < 40, is given by 

C(n,m) w c„n2™, (58) 

where 

ci = 0.03, C2 = 0.002, C3 = 0.0006, C4 = 1.5 x 10"^ and C5 = 3 x 10"*^. (59) 

There is a slight difference between even and odd cases, but this is insignificant in an estimate like this. The value 
of C5 is set to fit the observed trend in the other coefficients as the asymptotical behaviour is only glimpsed. In cases 
m = 2 and to = 4, it is vitally important to remember to truncate the asymptotical eigenvector ip^n^''^"^ correctly. 

For larger values of n, one needs very small values of tu in order to obtain accurate or even reasonable results. But 
for relatively small values of n, say n < 10, the error is extremely small at w w 0.01. The strong dependence on n 
means that the first few states can be obtained to a high precision even for quite strong couplings in the neighbourhood 
of a; w 0.1. We have determined the ground state n = up to the 31*^* order and numerical comparison strongly 
supports the asymptotical behaviour lu™ for m < 13. 

In order to make the above discussion more concrete, we explicitly give the second-order solutions as functions of 
uj, — ^/ujx, n (not n) and k' . For even values of n we find the solution function 

k' 

= Ar.^,, exp (-(i + (3 + 2n)uli2)e + (co/96)e^) ^ (/^i"^C'*^(l + (3fc - fc' + 10A:A:')w/48)) , (60) 

k=0 

where An^xo is a normalisation factor which ensures that || V'i^'^"'' II = 1- For odd values of n the result is nearly 
identical, i.e. 

k' 

^-(,"■2) = An,x, exp (-(i + (3 + 2n)u:/Z2)e + {uj/QQ)^) ^ + (7fc - fc' + 10fcA:')w/48)) . (61) 

fc=0 

The tiny difference Zk — > 7k in the generalised Hermite polynomial is very important, because otherwise the asymp- 
totical convergence 

||4'"^-o)_^(^.o)|| does not appear. The common exponential part in the third order solution 

function V'i"'^^ reads 

exp (-(i + (3 + 2n)tj/32 + (53 + 69n + 2W)uj'^ / IbiQ)^'^ + (w/96 + (11 + 6n)wVl024)C'* - (wVl280)C^) . (62) 

The explicit solution function f/'i"''"'* solves the asymptotical eigenvalue equation up to the order uj™" and yields a 
normwise convergence of ~ w™. 

When employing these asymptotical solutions, one should first study, how accurate eigenvectors are required for 
the problem at hand. The next step is to choose the order of the solution and the correct truncation with respect 
to X. Then, the calculations are performed and the results are obtained, hopefully faster than with the conventional 
approach of numerical diagonalisation^ 



IV. COMMENTS ON SOLVING THE ANSATZ 

In this section we discuss how to solve the set of algebraic equations resulting from Eq. (|33|l as effectively as 
possible. First we observe that the zeroth order, i.e. terms proportional to arbitrary powers of ^ are satisfied by the 
fact exp(O) — 1. Next, all equations related to terms 

{w^"+2-2''}i+fc' (63) 

are identically satisfied because of the recursion relation (|SJ) rewritten in terms of the coefficients h^jf^ . A careful reader 

notices that terms proportional to /3{,"^^2 appear, but they identically cancel and thus they are not constrained in 
this order. 
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From here on, we proceed by recursively solving the coefficients for the next order and also for sufficiently many 
values of n so that all coefficients in the expansions of {a^'^''} ^^'^ il^ki"^} have been constrained. In reality, we first 
obtained the solution function 'ip^~^'"^~^^ and a poorly formulated expression for arbitrary second-order solution, i.e. 
^(n,m-2) ^ but let US proceed in the way this should be done. Because the equation are quite difficult to handle with 
pen and paper, we chose to write and simplify the equations with Mathematica software.''^ 

We first consider the cases n = and n = 2 as simple examples. For n = we expand Eq. up to and including 
order oj'^ to find 



uj_ uj'^{x^ + 1/16) 
"2 ^ 2 



X 

~2 



1 

"48 



1:2 



ai,3 



3 1 



512 



Immediately, we obtain 



(0) 
"1.2 



-3/32 and oi'2\ 



1/96. 



-i-2. 



(0) 
1.2 



6a 



(64) 



(65) 



Inserting these into Eq. H64(l gives 



5 a 



(0) 
1,3 



-53/1536. 



In the case n — 2 and we examine all terms below the order of ui'^ . The generalised Hermite polynomial now reads 



Expanding all terms and moving them onto the same side yields the equation 

521 



(66) 




37a 



(2) 
2,2 



«S ( f - mfl 



(67) 



Notice that all terms proportional to llP and have canceled out, which again shows that the lowest-order approxi- 
mation for the eigenvalue and eigenstate are already correct and agree with the results for the continuous case. The 
three coefficients related to the -i/i^^^^) can be solved from the coefficients of w^, lS^x^ and lo'^x'^ and they read 



(2) 
"1,2 



-7/32, 



(2) 
"2,2 



1/96, and 



M'] = 1/4. 



(68) 



Substituting these into the set of equations and extending the calculation to order we find the subsequent coefficients 
to be 



1,3 



-275/1536, 



4^1 



= 23/1024, 



4^1 



= -1/1280, and 



37/256. 



(69) 



tj,; appear as constants, not unknowns. 



After solving a sufficient number of coefficients a^"-* and Z?^"-* one should start searching for regularities in the solution. 

Almost immediately we guessed the polynomial character of a^"'' , first in terms of n and later noticing that they 
should be written in terms of h as in Eq. H44() . This considerably helps solving the coefficients /S^"'' as for larger values 
of n the coefficients a^"^ 

In the beginning, we tried to solve all possible terms up to a given order in uj. First one should notice that only 
terms with I < m are required for the solution function 'ijj^x^'"^\ Assuming that the previous orders have been explicitly 
obtained, means that only the equations corresponding torn' = m in Eq. (|34|l have to be solved. In addition, generally 

in) 

known coefficients a). ^ identically satisfy equations corresponding to the highest powers of ^. Explicitly, if we assume 
that coefficients {"i' mlfe^fco ^-""^ known, only the equations for 



r m /:n+2ko — 2l' \kQ+k' 



(70) 
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are required and the expansion of Eq. (|33|l has to be carried out up to and including the order ^^^"i+'^o+n/a j-qj. 
coefficients I < m. 

After obtaining a rather complicated expression for the coefficients /^^"Lg, we happened to transform it into form 

equivalent to the present form and conjecture the general form of /3^"' in Eq. H51|) . The most important lesson taught 
by the discretised harmonic oscillator when solving the coefficients is that your numbers may be wrong, but the 
general forms usually are not. On several occasions, this became painfully obvious when the numbers did not check. 
Each and every time the general forms were correct, but the used expansion of Eq. H33() or the numbers inserted into 
it were not. 

Later on, we started to study the regularities in the general expressions. The polynomial structure of the coefficients 
a|.';' that do not contain any Gamma functions was relatively easy obtain, but the other set required a real stroke 

of luck. We managed to write some of these coefficients aj.^ as explicit products. After being pointed out, by 
Mathematica, that the first two could be written in terms of Gamma functions, it was only a question of finding the 
correct Gammas before Eq. H45() was written. In order to appreciate the technical part of obtaining the general form of 
the coefficients we point out that the coefficient oe'j!^l+4 was completed by solving the 12"^ order solution V'i™o ^^'^"^ 
confirmed by the case n = 1. Further terms have been obtained by solving the asymptotical differential equations H35|l . 

The regularities in the coefficients {/Sfc""*} have been found out using by studying the expansions with respect to k 
and k' . By conjecturing the recurring appearance of (lOfc' — k) in Eqs. H52|l and (|53() it becomes possible to solve the 

quantities defined in Eqs. H54 |l - H57() . In addition to these, the general expression for p1\_^ can be obtained from the 
known coefficients. 

Finally, we will estimate the difficulty of obtaining the explicit asymptotical solution We assume that both 

the expansion of the eigenvalues up to the required order and the solution ipi™ ^''^"^ have been obtained in advance. 
The coefficients a[,^ can be determined from the exponential parts of the eigenvectors up to and including the case 
n — m ~ 1. The completely general expressions in Eq. (|45|) are finished at much slower a pace. The asymptotically 
satisfied differential equations speed up this process considerably. 

Obtaining the coefficients /3^^ is more difficult. The general form H51|) shows that all states up to n = 4m — 3 must 
be solved. The explicit expressions for the leading parts, i.e. B{1, 21 — 2), B{1, 21 — 3), B{1, 21 — 3), and B{1, 21 — 4) 
make this task easier by 5 states. Thus all states up to n = 4Z — 8 must be found, unless further general properties 
are found. 

Regardless of these simplifications, the number of required terms and participating equations grows quite fast. 
Obviously, the general form of the coefficients in the exponential factor is much easier to obtain and thus they should 
be applied as early as possible. It is also possible that considerable simplifications or generalisations for the known 
coefficients lurk just around the corner. This has already happened on several occasions so far. We still choose to 
pause here, as the given general expressions have been validated rather convincingly and it not obvious, how, if at all, 
the next orders in the expansion would improve the results qualitatively. Wc hope a solid foundation has been laid 
for those striving towards the complete, asymptotical solution for the discretised harmonic oscillator. 

V. PROVING THE SOLUTION AND SOME GENERAL PROPERTIES 

Finally, we attack the difficult problem of actually showing that the solution is a general one. Thus far we have 
solved the equations for an increasing number of eigenstates using Eq. H33|l. This formulation is the best if actual 
numerical values of the coefficients a^,""* and are sought after. This is explained by symbolic math being most 
effective when the number of unknowns and symbols is as small as possible. In principle, the process explained below 
could be used for obtaining recursion relations between the coefficients of the solution and, subsequently, the full 
solution. Presently, we only show that the equations corresponding to leading orders up to are satisfied identically. 

We have to solve the equations corresponding to ^ijj">'^^+"^~'^i- in order to obtain the m}^ order solution. We 

have now obtained the explicit solution up to the seventh order so we can check if it is correct. For this purpose, we 
must write Eq. H33|) explicitly in terms of u := y/uj and i^, although odd powers of u eventually cancel. Multipliers of 
a^"^ and now read 

respectively. On the right-hand-side the non-trivial term is given by cj^^/2. Expanding all terms multiplying a fixed 
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term ft,^"^ up to the order yields terms 



4"^ 



2{k-k')-n/2 + 



^2 {n-l + 2{k-k')){n + 2{k-k')) 



C2 



= ftrMi + « + y ) )• (72) 



2 ^2 

The terms proportional to cancel and equating each power of ^ separately yields an equation 

h^^h{k' - fc) + [2{k + If T{k + l)]/ii"\ = 0, (73) 

where the signs + and — corresponds the even and odd values of n, respectively. The above equation is identically 
satisfied by the Hermite polynomials, which proves that the first order solution ipn'^"'^ is correct. A careful observer 
immediately asks about the second order corrections /3{,"2 which also yield terms proportional to . However, these 
coefficients are not fixed at all by Eq. (|33|l in the order v? . The only term that is easily solvable from this relation 
in the dominant coefficient = —1/2 which removes from the recursion relations. Later on, the dominant 

coefficients {cx^^}^^^ cancel the term h!"^}^^_^ in the order w™. 

In the next order we insert the solved coefficients and obtain for even values of a recursion relation 

6(n + 2 - 2fc)4"\ + [(2fc^ + fc^(42 - lln) - 6n - i-a^ - k{-6 + 9n- 5n^)]h["'^~ 

(1 + fc)(l + 2fc)(22 + 31k + fc2 - 5n - 5fcn)4+i + 2(2fc + 4)(2fc + 3)(2fc + 2)(2A: + 1)^+2 = 0' (74) 

which is again identically satisfied by the Hermite polynomials. For odd values of n, we find a similar recursion 
relation, once we replace k' — n/2 by k' ^ [n — l)/2. This completes the proof in order lo^ and validates the 

second-order eigenvectors ij^n'^"^ ■ 

In the third order the recursion relation for even values of n reads 

180(-4 + 2/c - n)/i^"'2 + 30(62 - A2k - 24fc2 + Ak^ + 74n - lOfcn - 22k^n + 1771^ + lO/cn^)^"^-!- 
[-450n - 450^2 - 90n^ - IQk^ + A:''(-452 + 105n) + ^^(-2332 + 2458n - 2,0Qn^)+ 
/c2(4230 + 4912n - 120471^ + 125n^) + A:(-300 - 585n - 1258n2 + 179n3)]4"^ + 

(1 + fc)(l + 2fc)(1022 + 2705A: + 3684fc2 + 326/fc^ + 5fc* - 2274n - 4430fcn - n2Qk^n - 50k^n + 45471^+ 

579kn^ + 125fcV)4"\ - 16(1 + k){2 + k){l + 2fc)(3 + 2A:)(110 + lOl/c + 5k^ - 50n - 25A:n)4^2+ 

256(1 + k){2 + fc)(3 + fc)(l + 2/c)(3 + 2fc)(5 + 2k)h'i!'^^ = 0. (75) 

Because the Hermite polynomials satisfy this and the corresponding relation for odd values of n the solution ipn'^""^ 
has been rigorously proven as correct. 

The eigenvectors ^/ji™'^"-* tend to the eigenvectors ipl^'"^ of H(xo) at an asymptotical rate proportional to lu"^. The 
exact eigenvectors are orthogonal as eigenvectors of a Hermitian matrix and by their closure relation we can write 

^^-,.0) ^ ^(.0) + ^rn J2 6„,vi?\ (76) 
n' 

where 6„' are finite constants such that ^„ l^nP < oo in the limit lj — )■ 0. The orthonormality relation for the 
asymptotical solutions thus reads 

(V'i'"'^°)|^i"'^"^)=5w+0(c^'"), (77) 

provided that the sum |6„| is finite for both states. In other words, the eigenvectors T/ii"*'^"^ become orthonormal 
at the rate of w™. Numerical checks seem to confirm this, at least for relatively small values of n. 

As a final effort, we outline a plausible "proof for the asymptotical convergence. As a first step, we show that with- 
out loss of generality we can examine a finite truncation of the eigenvector i/'i"'^"'', the vector V'i^jo •= {V'j-'Lo J'j=-jo 
for sufficiently large jo- For sufficiently large values of |j| the eigenvalue A„ becomes insignificant in Eq. H32|) and we 
write an approximate equation 

(V'i"_\-24")+^i;\)/(x2)=c.2. (78) 

For sufficiently large values of x and/or j the sign of ipx^^ is constant and this equation shows that the function 
-ipx = exp(— ^(a;^ ~ Jo)/(2 + £))j for some small e > 0, is a dominant sequence for Now, the limiting 

sequence of norms 



lim (79) 



71, IQ 
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vanishes exponentially with respect to Jq. In other words, we can always find a finite jo such that the error in the 
norm is sufficiently small. 

Next we use the fact that solution satisfies the eigenvalue equation (p!^ up to the order w™ when written 

in terms of ^. Thus we can write 



( (n,m) , ( (n.m) 

7"",^ "^"^^ = l + 0(c.™+i). (80) 

2#''"^(-A„+c.e/2) 

We fix the scales of the eigenvectors by setting (V'rr''^°^)j — {'>Pri°'')j for an arbitrary j. It would be very tempting 
to say that Eq. (|80(l implies — ('(/'i^"'')j+i(l + 0{uj™^^)) and then wonder why convergence is not 

asymptotically proportional to lu"^^^. As already explained the solution ■0^™'^°-' does not fix the coefficients /^["^j^^ 
which most definitely yield terms proportional to w™. Thus we obtain a relation 

= (V'^^Will + Oicu"^)). (81) 

By matching the eigenvectors at j = 0, and expanding the components to the finite values ±jo shows that the order 
of error is w™ for all components with \j\ < jg. Becase the error caused by the truncation is insignificant the result 
holds for the full eigenvectors and we obtain the desired result 

|,^(™,.o)_^(xo)|, (82) 
or at least show that the result is quite plausible. 



VI. CONCLUSIONS 



We have obtained an explicit, asymptotical solution for the discretised harmonic oscillator. Both the eigenvalues 
and eigenvectors have been obtained and we can choose a prespecified rate of convergence towards the exact solutions. 
This is done by truncating the ansatz solution accordingly. Because the problem can be mapped onto the Mathicu 
differential equation, we simultaneously provide asymptotical expressions for the Mathieu functions. The Schrodinger 
equation of the quantum pendulum corresponds to the Mathieu equation, which yields immediate applications for the 
results. 

The method described above can be generalised to accommodate several coordinate dimensions with only minor 
changes. This should make the results of Ref. 29 both more transparent and more rigorous. The tunncUing-charging 
Hamiltonian of a Cooper pair pump corresponds to a modified multi-dimensional Mathieu equation. 

Alternatively, ansatzes similar to Eq. H30|l could be constructed in case of difference equations that become identical 
to analytically solvable differential equations in some asymptotical limit. Initially, the problem assumes the form 
of an infinite-dimensional, two-parameter [eigenvalue] problem, where the asymptotical solutions [eigenvalues and 
eigenvectors] must be obtained. The ansatz maps the problem onto an infinite set of algebraic equations that must 
solved. If the form of the ansatz is correct, one may determine some general properties of the exact solution. 
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